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Abstract 

The D°-D° mixing is analyzed in a weak gauged U(4)l <8> U (4)^ chiral lagrangian 
model where the electroweak interaction is introduced as a gauge theory over the 
meson degrees of freedom. This model allows a particular realization of the G.I.M. 
mechanism and then could be useful in the study of processes where G.I.M. suppres- 
sion is effective. As a test of the model we have also analyzed the K°—K° mixing. 
We find ArriK in good agreement with the experimental result and we show that 
the D°-D° mixing is very much suppressed in agreement with previous estimates 
in the Heavy Quark expansion framework. 



1 Introduction 



Higher order effects in the perturbative treatment of the electroweak Standard Model 
(SM) play an active role both in testing some of its fundamental ingredients and in 
uncovering phenomena not explained by the Standard framework. One of the peculiarities 
of the SM is the absence (at leading order in the perturbative expansion) of flavour 
changing neutral currents (FCNC). This characteristic feature strongly affects flavour 
mixing and rare meson decays. 

Neglecting the tiny CP-violating amplitudes (as we do in this Letter) the perturbative 
contributions to P°-P° mixing are straightforward to evaluate in the SM. However the 
numerical quantification is affected by our rather poor knowledge of some of the relevant 
elements of the CKM matrix and of the quark masses. In any case more dubious is the 
treatment of long-distance contributions. These have two main sources: a) the evaluation 
of the matrix element of the four-quark AF = 2 operator in the effective hamiltonian (F is 
short for flavour), and b) the handling of mesonic intermediate states contributing to the 
dispersive amplitude. Both features are characterized by a manifestly non-perturbative 
origin. 

The relative importance of short vs. long-distance contributions depends crucially of 
the flavour involved. This is so because as the flavour is heavier the transitions involve 
higher transfer of momenta and therefore non-perturbative corrections are less relevant. 
As a consequence B®-B® and B®-Bj are expected to be dominated by the perturbative 
regime. This is not so clear for K°-K° and definitely not the case of D°-D° (due to the 
possible contribution of resonances in the highly populated kinematical region of charmed 
mesons) . 

In this Letter we address the issue of non-perturbative corrections to the mixing and 
therefore we will be concerned with these last two cases. 

At hadronic level meson mixing with AF = 2 transitions has ushered permanent inter- 
est because its implications : tests of FCNC, close relation with CP-violation, prospects 
of new physics beyond the SM, etc. (see [|l], ||. |3|. |J and references therein). Mixing occurs 
through radiative corrections in the SM. The mass eigenstates are 



and the mixing produces P° «-> P° oscillations of amplitude exp(iAMt) with AM = 
m 2 — rrii — 2(r 2 — T 1 )/2. The parameter controlling the oscillation is x = Am/T [] and 
consequently the mass difference between the eigenstates. With our normalization this 
mass difference is related with the AF = 2 transition through 




(2) 



1 We call x k = Amx/rs and xd = Arra^/rc. 
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where mp is short for the relevant pseudoscalar mass. 

K°-K° and D°-D° mixings have notable different features in the SM. Both mixings 
have a complementary interest because test FCNC in the upper and lower sector of the 
families, respectively. Experimentally f| xk — 0.5 and \xd\ < 0.09 and, as we will shortly 
see, the perturbative evaluation in the SM is able to point out these different behaviours 
because it predicts Am K ~ m c and Am D ~ m s . Moreover the perturbative evaluation 
gives the right order for the measured K°-K° mixing while it gives a very small D°-D° 
mixing. This is the case because GIM mechanism || is much more effective in the lower 
sector of the families than in the upper one due to the quark mass differences. 

Dispersive amplitudes due to intermediate mesonic states could give relevant non- 
perturbative contributions. These have been considered by several authors 0, [7], PL [II], 
|TTf and the conclusion is that they indeed are important. In D°-D° could even give a 
result several orders of magnitude bigger than the perturbative one. This might be due to 
the fact that the expected GIM suppression could be spoiled by the bad SU (3) breaking 
observed in several decays of charmed mesons [[], . 



We have proposed a weak-gauged chiral model up to and including charmed mesons in 
order to implement at hadronic level all the structure of symmetries of the SM, specifically 



including the GIM mechanism JT3|. Therefore, the SU(3)l <8> SU(3)r chiral symmetry 
of the strong interactions must be extended to a U(A)l ® £/(4)# symmetry in order to 
implement the weak SU(2)l representations. Nevertheless it must be clear that our model 
is not pretended to be a model for the strong interaction. Of course chiral symmetry is 
badly broken in the charm sector. We expect that this model may provide a realistic 
approach to those processes where relevant features of the flavour sector of the SM like 
GIM play a significant role. In this Letter we evaluate in our model [14|, the leading order 
result for the K°-K° and D°-D° mixings. 

In Section 2 we remind briefly the perturbative evaluations of the mixings in the SM. 
An outline of the basic ingredients of our model is sketched in Section 3. Then in Sections 
4 and 5 we specify our results for K°-K° and D°-D° , respectively. We end with our 
conclusions in Section 6. 



2 The perturbative regime 

The short-distance contribution to meson mixing in the SM is represented by the 
usual box diagrams. In order to carry out this calculation one starts with the weak 
effective hamiltonian 'H^ff =2 in terms of four-quark operators. We explain both K°-K° 
and D°-D° mixings in turn. 
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2.1 K°-K° 



The effective hamiltonian (we are not interested in CP-violation) including the next- 
to-leading logarithmic QCD corrections, is given by |L5|, |l6f 



n_/AS=2 
H eff 



Gp ml X 2 C r]i S 



16tt 



a 



(/i 2 )( -2/9) Q AS=2 



where 



O 



AS=2 



A c = VcdVcs with Vij the CKM matrix, a and f3 are colour indices and 



S{x) 



1 + 



9 
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Qx 2 



x 



1-x) 



1 



lnx 



x) 



(3) 



(4) 



(5) 



In Eq. (|^) 7]i carries the information of strong QCD corrections flT7| . 

A part of the perturbative uncertainties in order to get a firm prediction we need to 
do the evaluation of the matrix element of the four-quark operator in Eq. @ between the 
asymptotic kaon states. This, of course, involves the unsolved problem of hadronization 
and our ignorance is expressed by the B parameter defined as 

' ! ' " (6) 



(W\O AS = 2 \K ) = -fW K B{v) 
where, in the leading logarithmic approximation in the strong coupling, 

B(Ai) = B K a s {^) +2 " . 



(7) 



A simple factorization approach and vacuum insertion |T| gives B{^i) = 1. The non- 
perturbative parameter Bk has been calculated using lattice gauge theories 1/N C 
expansion and QCD sum rules 2(J leading to Bx = 0.7 ± 0.2. The QCD hadron 
duality approach |2l| favours lower values B K = 0.4 ± 0.1. 



2.2 D°-D° 

The effective hamiltonian (neglecting the contribution of the b quark due to the sup- 
pression of the relevant CKM matrix elements) is given by |3], |9| 

H%f 2 = --^G 2 F \ s \ d K ~f df ( 0* c = 2 + 2 0* c = 2 ) , (8) 
where \ = V*Jf U i and the four-quark operators are defined as 

(9) 

o$ c = 2 = (bZc£)(b£c£) . 
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This last operator has no analogous in the K°-K° case and it is due to the non-negligible 
charm quark mass m c carried by two of the four external legs of the box diagram. 

Analogously to the K°-K° mixing the evaluation of the matrix elements of the two 
four-quark operators involves non-perturbative aspects that, this time, are parameterized 
by two quantities Br, and B' D for the two operators in Eq. (^) respectively, and defined to 
be the unity if factorization and vacuum insertion approximations are used. Corrections 
to this value are potentially large but are not expected to change the order of magnitude. 

At any rate from Eq. (|8|) we can conclude that D°-D° mixing is suppressed typically 
by a factor ~ m\jm\ over the K°-K° case. 



3 The model 



Our lagrangian is [13| : 

£ = £ 

mesons + £-Higgs + £hM + • • • , (10) 

where the dots are short for pure gauge boson terms. Here £ me sons contains the strong 
interaction between mesons and the couplings of mesons to the gauge bosons. We have a 
set of 16 scalar and 16 pseudoscalar fields which we assign to the (4, 4) © (4, 4) represen- 
tation of the chiral {7(4)/, <g> U(4) R group. We denote the meson matrix by U = £ + ill, 
where £ is the scalar and II the pseudoscalar matrices of fields. The explicit expression 
for the pseudoscalar matrix is 



' 2 r l° 1 ^2 



n 



75% + ^15 



7T 



K- 



71 



In _ J_ 71-0+ 

2'l° V2 ^ 



75% + 7^15 



K° 



K° 



ho - 76%+ 



D° 



D~ 



d: 



12 %5 

D 



A similar matrix can be written for the scalar mesons. Our notation for scalar mesons 
is cr ,a 8 , <j 15 , cr + , o" 3 , /t, 6, 5s instead of i] ,T] 8 , ^ 15 , 7r + , 7r°, D and D s respectively. 
With these definitions £ meS ons is 



Imesons = -Tr[{D^U')\D a U')} - V chlral (U) 



(12) 
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where V c hirai is the chiral potential 



VchiraliU) = -fM 2 Tr(U*U) + 

H 2 [a Tr{U ] U) 2 + b (Tr(U ] U)) 2 + c (detU + detU ] )} , 



(13) 



with fil > in order to develop spontaneous breaking of chiral symmetry. The covariant 
derivative is: 



DfJJ' = 8,U' - igf ■ %U' - ig'YiB^U' + tg'U'Y R B^ 



with 



U' = SUS ] . 



(14) 



(15) 



Here and are the gauge bosons related with the SU (2) l and the U(l)y groups. The 
7^ matrices are the SU (2) generators and Yl and Yr are the left and right hypercharges. 
The matrix S will be the Cabibbo rotation once the SU(2)l <8> U(1)y symmetry gets 
spontaneously broken. With our definitions for U we have 



f 

T X TT X 



Yr 



4x4 



Y, 
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\ 



cos^c sin 9 C 
— sin 6 C cos 9 C 
V 1 / 



(16) 



with t the usual Pauli matrices and 9 C the Cabibbo angle. The charge operator is Q = 
Y R = Y L + T 3 . 

In Eq. (0) Cm ggs is the usual lagrangian for the minimal model of Higgs of the 
Standard Theory. Chm, finally, is a Higgs-mesons coupling term which will give masses 
to the mesons after the spontaneous symmetry breaking of the weak symmetry. In order 
to write in a compact form this term we introduce the usual Higgs doublet as a 4 x 4 
matrix: 

I -fckl>-ix) o 



H 



V 



— s~ 





s + 






-s~ 



\ 



(17) 



We consider only the simplest local gauge invariant term 



C HM = Tr(AS ] H ] SU + h.c. 
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where S is given in Eq. (16) and the most general form of A assuming isospin symmetry 
is 



.4 



V 



a 



a 



7 



\ 



8 ) 



(19) 



Chm breaks explicitly the S77(4) ® S77(4) chiral symmetry. Through the spontaneous 
breaking of the chiral and weak symmetry the matrices of mesons and Higgs get a non- 
zero vacuum expectation value 



(°\U\o) = F = 



( u 



V 



./;, 



fs ) 



(o|if|o) = -l=0 o J 4x4 



(20) 



Therefore from Eq. (EH3) we get 



HM 



^j=Tr{A{U + U ] )) + Tr{AS ] H ] SU + h.c] 



(21) 



where 



H = H 



(o\H\o) . (22) 

We note that the lagrangian in Eq. fll8|) transforms under chiral SU(4) <S> SU(A) as the 
(4, 4) © (4, 4) representation and, therefore, the first term in Eq. ( |2"T| ) is similar to the 
usual explicit breaking of chiral symmetry. 

Some aspects of this model are worth to emphasize: 

a) The starting point of our ideas is to believe that the symmetries of the Standard 
Theory are essential to describe the weak processes of hadrons, and any model for 
them has to support not only its symmetries but its symmetry breaking patterns 
too. Our realization with two complete families allows to implement rigorously the 
structure of the SM at hadronic level. 

b) The GIM mechanism is naturally implemented in our scheme. Therefore we have 
not flavour changing neutral currents at leading order in the perturbative expansion 
and in this way the model might presumably be trusted in the study of processes 
where GIM plays a significant role. 

c) The model has been tested satisfactorily in the study of non-leptonic D decays into 
two pseudoscalars [22], non-resonant D — > PPP and also, at one-loop level, 
in radiative rare kaon decays (K + — > n + £ + £~ 
results. 



K, 



it 



||13| with satisfactory 
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Hence we think it is worth to use the model just described in order to study K°-K° and 
D°-D° mixings. We have already seen in Section 2 that the short-distance contribution 
to K°-K° mixing is dominated by the charm degree of freedom while the strange quark 
dominates D°-D° mixing. Our model provides these ingredients and, moreover, allows the 
evaluation of long-distance dispersive contributions mediated by the pseudoscalar (and 
maybe scalar) mesons. 

Therefore we proceed to calculate them at 0(G 2 F ) in the weak perturbation expansion. 
Due to our ignorance on the scalar masses and in order to present simplified analytical 
expressions, we first analyze the case in which m sca i ar 3> m pseu doscaiar, that is /j% — > oo in 
Eq. ( |13|) but with /^c = d constant in order to keep the r) mass finite. In this case F in 
Eq. (EDI) is 



+ O 



m 



Pseudoscalar 



4x4 T ^ 2 



(23) 



with 

f ° = ^hi ■ < 24) 

With our normalization / = f n ~ 90MeV in the chiral limit. 

Using this expansion, at leading order, we get a relation between the pseudoscalar 
masses, 

2 2 2 2 / n r\ 

m D s - m D = m K - m n . (25) 

Eq. fl2"3"|) is a crude approximation for the charm sector, we must remember here that 
U(4:)l <8> U(4)r is not a symmetry for the strong interaction in nature. Therefore we have 
also considered the case of finite scalar masses in both calculations where one goes beyond 
the rough approximation represented by Eq. (^3|) and ([25]) . 

We use the renormalizable gauge |24| where there is no mixing meson-W and 
instead we have a meson-charged Higgs and meson mixings as shown in Fig. 1. 

At this order the two different topologies of Feynman diagram contributions to the 
mixing are those in Fig. 1. Both diagrams are separately finite but only the sum of the 
two is gauge independent. 



4 K°—K° mixing 

In order to test the model we can calculate Am^. The problem of evaluating long- 
distance contributions to the K°-K° mixing has been addressed extensively in the bibli- 
ography [0, 10, 25, p6[] . These are due to intermediate hadronic states in the oscillation 



and are presumably dominated by the octet of pseudoscalars: K° <-> (7r°,r]) <-> K°, 
K° «-> 7T7T <-> K°, etc. An estimate of these contributions gives results comparable with 
the experimental mass difference. However the evaluation has significant uncertainties. 
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Figure 1: The loop diagrams are the contribution of the model to the K°-K° and D°-D° 
mixing. The crossed propagator is the sum of the Higgs s contribution plus the meson 
mixing. 



The computation of the K°-K° mixing in our model through the Feynman diagrams 
of Fig. 1 gives for the difference of masses: 
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"I 2 2\2n/ 2 2\ , / 2 2\( 2 2\r>t 2 2\ 

-jK - m J p {m^ m J + ( m K ~ m D)i m K ~ rnJ p ( m D> m J 



where P(m 2 ,m2) is defined through 



I 2 (m 2 ,ml) 



167T 



{ + 7 - ln(47r) - 2 + P(m 2 ,m 2 ) } , 



(26) 



(27) 



and l2{rn\,m%) the scalar two-point function. 

Several comments are in order about our result Eq. (|26|): 

- Our result is finite and only depends on the masses of the pseudoscalar mesons and 
the meson decay constants once the mass expansion induced by the assumption 

TTlscalar ^ m pseudoscalar IS Carried Out. 

- We can easily see the effect of the GIM cancellation (in the SU(4)f limit m u = m c ) 
by inputting mu = m n and mx = ^d s in Eq. that gives evidently Amx = 0. 
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- From the observation that the short-distance contribution to Arrix scales with m 2 
in Eq. (Q) one should expect in our model that our result for Am^ might scale like 
m\). However using the mass relation Eq. ( f25|) in Eq. ( p6|) we see that the term in 
mjj cancels. This fact provides a further suppression beyond GIM mechanism. 

Using the mass relation in Eq. (EST) and expanding the logarithms in turn the result 
in Eq. (|26| ) can be written as 

1 1 9 

Am K = —^G 2 F sin 2 9 cos 2 9 fl — (1 + u K ) m 2 D m K . (28) 

47T 11 

(m 2 \ 

where uj k = OR5 and our result gives ~ —0.09 providing a tiny 10% correction 



ni 



D 



to the leading term. Using the values of the decay constants and masses of || we get 

Am K = 3.45 x 10~ 15 GeV . (29) 

Inputting the experimental result for T s @ we predict xk = 0.49. Our results are in 
rather good agreement with the experimental figure || 

Am" p = (3.491 ± 0.009) x 10~ 15 GeV . (30) 

This fact can be understood because the dispersive part of the amplitude is dominated 
by the AI = 3/2 transitions and, therefore, through the (7i + 7T \T-lw\K + ) amplitude [p5 
that is well recovered in our model. 

One could also consider the corrections coming from finite scalar masses. These happen 
to depend only on two parameters : / 7 // Q defined in Eq. (|^) and d related to the 
masses of the mesons rj and v( . This dependence is showed in Fig. 2. Inside reasonable 
values of these parameters our prediction Amx in Eq. (|29| ) could only be modified in a 
factor 2. We can conclude that the model gives reasonable values for Arrix, in the range 
(1.8-3.5) x 10- 15 GeV. 



5 D°—D° mixing 

As has been noted previously in Section 2 the short-distance contribution to the D°-D° 
mixing is very much suppressed relatively to the K°-K° mixing, roughly a factor m A s /m A c . 
A careful evaluation shows that the suppression is not that much but in any case the 
perturbative contribution gives 

Am^ ~ 10- 2 Am^ x . (31) 

The study of long-distance contributions to D°-D° mixing has been considered previously 
|| ||, [Tl[| with very different conclusions. While Wolfenstein , Donoghue et al. || and 
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Kaeding [12] conclude that the dispersive non-perturbative contribution to D°-D° mixing 



might be huge at the level of providing Amo ~ Am^ due to the fact that there is a 



strong SU(3) breaking that overcomes the GIM suppression, Georgi JTT] using arguments 
of a Heavy Quark expansion concludes that this is unlikely, and a careful study in this 
framework |27) seems to confirm this assessment. 

The evaluation of Amp at leading order in our model goes again through the corre- 
sponding Feynman diagrams in Fig. 1. By expanding in m Pseudoscalar /m Scalar and keeping 
the leading term we get 



G F sin cos Of r 2 ■_> , 2 , 2 2 , , 2 1 2 2 1 2 > 

Am D = —5 (m K - mj + [m K - mJ(m D In m D - m D \nm D ) 

S7r rrif) 



1 2/2 1 2 



1 2 / 2 1 2 
+ (m n + m D 




1/ 2 

~o( m ^ ~ m l) 2 p ( m l^ m l) + ( m D ~ m K)i m D P(m 2 K ,m*) 

(32) 

where P(mf,rn|) has been defined in Eq. (p7|). Again we find a finite result in terms 
of the pseudoscalar masses and the meson decay constants. Moreover GIM suppression 
(for exact SU(3) P symmetry, m d = m s ) is translated now into the meson language as 
rriK = rn-K and mc s = mp that gives Amu = 0. 

Using the mass relation in Eq. and expanding over the charmed meson mass we 
have: 

Am D = -^G% sin 2 cos 2 fl f 1 - — (1 + cu D ) , (33) 

with cjd ~ 0.05 carrying the contribution of higher order in m\^jm D . As we see the 
Amp, scales with the mass of the kaon to the fourth over the mass of the charmed meson. 
We notice that our result implies a suppression 



Amp ^Am K . (34) 

m D 

From Eq. ([33|) we predict 

Amp = 2.21 x 10- 17 GeV , (35) 

and with the experimental |j] we have xd — 1-4 x 10~ 5 . There is only an experi- 
mental upper limit on Amp || 

\Am e * p \ < 1.38 x 10- 13 GeV (90%CL) , (36) 
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Figure 2: x — Amp/T as function of r = / 7 // a - This ratio is related to the mass of the 
scalar mesons. The left-hand side (r = 1) of the plots correspond to the limit m sca i ars — > 
oo and the right-hand side to m sca i ars ~ 2 GeV. The parameter d = —15, —30, —45 is 
related to the mass of the rj' meson; we indicate the corresponding mass in brackets. 
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that is four orders of magnitude bigger than our result. 

However, there are many dynamical effects in the energy region of 1 — 2 GeV which 
could influence the final states in D decays and then modify our prediction. Analogously 
to the K°-K° mixing we can give an estimate of the dependence of Am^ on the masses 
of the scalar mesons (Fig. 2). In this case the dependence is bigger than in Amx but for 
reasonable values of / 7 / f a and d there is no change in the order of magnitude estimate. 
Moreover we notice that the sign of Amp could be changed due to the contribution 
of the scalar mesons. We think that, in the SM, a bigger order of magnitude for the 
mixing than our result implied by Eq. (|35|) could only be expected through a 
large SU(3)p breaking possibly induced by final-state interactions generated through 
vector meson resonances. This effect we know it is important in non-leptonic decays of 
charmed mesons |28| and therefore could be relevant in the D°-D° mixing. However, 
it can be expected a cancellation in the total contribution of the vector mesons due 
to the GIM mechanism. Actually we expect, from these new intermediate channels, 
contributions without changing noticeably the order of magnitude. Once incorporated all 
the contributions, the final-state interactions correspond to a unitarity transformation 
and then it is not expected a modification in the mixing. 

A result for Amp as given by our prediction in Eq. (^), if true, is out of the ex- 
perimental reach in the foreseen future. However an experimental result at level of 
Amu > 1CT 16 GeV could be a signal of new effects outside of the SM. An overview 



of the experimental techniques and problems has been given in |[29|| . The recent experi- 
mental determination of the doubly Cabibbo suppressed channel D° — > K + n~ ||, even 
still with a big error to reduce represents a good step forward in this direction. 



6 Conclusion 

We have studied the K°-K° and D°-D° mixings in a weak gauged U(4)l <8> U(4)r 
chiral lagrangian model that incorporates the GIM mechanism at hadronic level. We 
have got a result for Amx in rather good agreement with the experimental result and 
a prediction for Amo that shows a strong GIM cancellation and provides xd in the 
range (4, —1) x 10~ 5 . Our result is in good agreement with previous estimates in the 



framework of a Heavy Quark expansion [O, Other calculations giving largest Arrijj 



are not really in contradiction with our result. Reevaluating the estimate in |J with the 
actual experimental values for the widths we obtain more cancellation than the originally 
assumed in [[J. Using theoretical results for the widths from [22] and Eq. ([32]) we can 



also observe that the phenomenological expression used in |§ must be corrected, giving 

2 In (22) we predicted, with the model used here, this width and our result T(D° — > K + n~) = 3.3 x 
10~ le GeV is in very good agreement with the experimental figure T(D° — > K + n~) = (4.6 ± 2.2) x 
lQ- 16 GeV. 
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additional cancellations. From the present analysis and previous works on the subject we 
can safely conclude that a value of \xr> \ over 10~ 4 would be a clear signal of new physics 
in the charm sector. 
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